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ABSTRACT 

We  present  the  first  coordinated  Very  Long  Baseline  Array  (VLB A)  /  Very  Large  Telescope  Interferometer 
( VLTI  )  measurements  of  the  stellar  diameter  and  circumstellar  atmosphere  of  a  Mira  variable  star.  Observations  of 
the  v  =  1,  J  =  1-0  (43.1  GHz)  and  v  =  2, J  =  1-0  (42.8  GHz)  SiO  maser  emission  toward  the  Mira  variable  S  Ori 
were  conducted  using  the  VLBA.  Coordinated  near-infrared  A'- band  measurements  of  the  stellar  diameter  were 
performed  using  VLTI-VINCI  closely  spaced  in  time  to  the  VLBA  observations.  Analysis  of  the  SiO  maser  data 
recorded  at  a  visual  variability  phase  0.73  shows  the  average  distance  of  the  masers  from  the  center  of  the 
distribution  to  be  9.4  mas  for  the  v  =  1  masers  and  8.8  mas  for  the  v  =  2  masers.  The  velocity  structure  of  the  SiO 
masers  appears  to  be  random,  with  no  significant  indication  of  global  expansion/infall  or  rotation.  The  determined 
near-infrared,  K-band,  uniform  disk(  UD)  diameters  decreased  from  ~10. 5  mas  at  phase  0.80  to  M0.2  mas  at  phase 
0.95.  For  the  epoch  of  our  VLBA  measurements,  an  extrapolated  UD  diameter  of  0^D  =  10.8  ±  0.3  mas  was 
obtained,  corresponding  to  a  linear  radius  of  =  2.3  ±  0.5  AU  or  =  490  ±115  R0.  Our  coordinated 
VLBA/VLTI  measurements  show  that  the  masers  lie  relatively  close  to  the  stellar  photosphere  at  a  distance  of  ~2 
photospheric  radii,  consistent  with  model  estimates.  This  result  is  virtually  free  of  the  usual  uncertainty  inherent  in 
comparing  observations  of  variable  stars  widely  separated  in  time  and  stellar  phase. 

Subject  headings:  masers  —  stars:  AGB  and  post-AGB  —  stars:  atmospheres  —  stars:  late-type  — 
stars:  mass  loss  —  techniques:  interferometric 


1.  INTRODUCTION 

The  evolution  of  cool  luminous  stars,  including  Mira  vari¬ 
ables,  is  accompanied  by  significant  mass  loss  to  the  circum¬ 
stellar  envelope  (CSE)  with  mass-loss  rates  of  up  to  10~4 
yr  1  (e.g.,  Jura  &  Kleinmann  1990).  There  are  currently  a 
number  of  tools  that  can  be  used  to  study  the  stellar  surface  and 
the  CSE  at  various  wavelengths.  Optical/near-infrared  long- 
baseline  interferometry  has  provided  information  regarding  the 
stellar  diameter,  effective  temperature,  center-to-limb  intensity 
variation,  and  the  dependence  of  these  parameters  on  wave¬ 
length  and  variability  phase  for  a  number  of  Mira  variables 
(e.g.,  Haniff  et  al.  1995;  van  Belle  et  al.  1996;  Perrin  et  al.  1999; 
Young  et  al.  2000;  Thompson  et  al.  2002a,  2002b;  Woodruff 
et  al.  2004).  The  structure  and  dynamics  of  the  CSE  of  Mira 
variables  and  supergiants  have  been  investigated  by  mapping 
the  SiO  maser  emission  at  typical  distances  of  2-4  stellar 
radii  toward  these  stars  using  very  long  baseline  interferom¬ 
etry  ( VLBI)  at  radio  wavelengths  (e.g.,  Boboltz  et  al.  1997; 
Kemball  &  Diamond  1997;  Boboltz  &  Marvel  2000;  Hollis 
et  al.  2001;  Diamond  &  Kemball  2003).  Dust  shells  outward 
of  the  SiO  maser  emission  zone  have  been  studied  using  mid- 
infrared  interferometry  (e.g.,  Danchi  et  al.  1994;  Greenhill 
et  al.  1995;  Townes  2003).  Near-  and  mid-infrared  interferom¬ 
etry  have  also  been  used  to  study  the  warm  extended  atmo¬ 
sphere  of  Mira  variables  containing  gaseous  H20,  SiO,  and 
CO  molecules  (e.g.,  Perrin  et  al.  1999;  Mennesson  et  al.  2002; 
Ohnaka  et  al.  2005). 

Theoretical  models  describe  and  predict  photospheric  center- 
to-limb  intensity  variations  including  the  effect  of  close  mo¬ 
lecular  layers  (e.g.,  Bessel  et  al.  1996;  Hofmann  et  al.  1998;  Tej 


et  al.  2003;  Ireland  et  al.  2004),  the  location  and  dynamics  of 
SiO  maser  emission  zones  (Humphreys  et  al.  1996,  2002),  as 
well  as  gas  and  dust  components  of  the  stellar  winds  (e.g., 
Winters  et  al.  2000;  Hoefner  et  al.  2003). 

Multiwavelength  studies  of  the  stellar  surface  and  the  CSE 
using  a  combination  of  the  above  techniques  are  well  suited  to 
study  the  structure  of  the  CSE,  the  mass-loss  process,  and  the 
evolution  of  these  stars.  It  is  common  for  infrared  observations 
of  the  stellar  photosphere  or  the  circumstellar  dust  to  be  com¬ 
pared  to  radio  observations  of  circumstellar  masers  from  the 
literature,  and  vice  versa.  For  example,  Danchi  et  al.  (1994) 
compared  the  inner  radii  and  extents  of  dust  shells  as  deter¬ 
mined  from  Infrared  Spatial  Interferometer  (ISI)  observations 
with  estimates  of  the  photospheric  radii  and  measurements  of 
the  extents  of  the  SiO,  H20,  and  OH  maser  shells.  Similarly, 
Monnier  et  al.  (2004)  compared  high-resolution  data  on  the 
stellar  diameters  and  dust  shells  of  two  stars  ( VX  Sgr  and  NML 
Cyg)  from  Keck  aperture  masking  and  Infrared  Optical  Tele¬ 
scope  Array  ( IOTA)  interferometry  with  previously  published 
maps  of  the  SiO  maser  emission  toward  these  stars.  Such  com¬ 
parisons,  however,  are  somewhat  limited  because  of  the  inher¬ 
ent  variability  in  the  CSE  and  the  star  itself.  The  SiO  masers, 
for  example,  have  been  shown  to  have  proper  motions  and  in¬ 
tensity  variations  on  timescales  of  a  few  weeks  ( Boboltz  et  al. 
1997;  Diamond  &  Kemball  2003).  Similarly,  interferometric 
measurements  of  the  stellar  parameters  vary  with  stellar  vari¬ 
ability  phase,  as  discussed  above.  This  inherent  variability  of 
Mira  stars  and  their  surroundings  necessitates  a  more  coordi¬ 
nated  approach  to  multiwavelength  studies  of  these  objects. 

The  first  such  coordinated  multiwavelength  study  of  a  late- 
type  variable  was  performed  by  Greenhill  et  al.  (1995)  for  the 
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star  VX  Sgr.  Coordinated  measurements  included  both  mid- 
infrared  observations  to  determine  the  extent  of  the  circumstel- 
lar  dust  shell  and  VLBI  mapping  of  the  SiO  maser  zone.  More 
recently,  Cotton  et  al.  (2004)  compared  Very  Long  Baseline 
Array  (VLB A)  images  of  the  SiO  maser  emission  of  several 
Mira  stars  with  near-infrared  diameters  obtained  at  the  IOTA 
interferometer  (Mennesson  et  al.  2002)  and  with  literature  val¬ 
ues  of  the  mid-infrared  dust  shell  extensions. 

Here,  we  present  the  first  results  from  a  series  of  coordinated 
VLTI/VLBA  observations  of  late-type  variable  stars,  in  par¬ 
ticular,  for  the  Mira  variable  S  Ori.  Our  results  include  the  first 
radio  interferometric  images  of  the  SiO  masers  (VLB A)  and 
the  most  accurate  photospheric  diameter  measurements  using 
near-infrared  A'-band  interferometry  (VLTI-VINCI). 

2.  CHARACTERISTICS  OF  S  ORI 

The  target  of  our  observations,  S  Ori,  is  a  Mira  variable  star 
with  spectral  type  M6.5-M9.5e  (Kholopov  et  al.  1985)  and  V 
magnitude  8.4-13.3  (Pojmanski  2002).  Whitelock  et  al.  (2000) 
report  a  quite  low  near-infrared  ./-band  pulsation  amplitude  of 
0.55  mag.  They  also  derive  a  mean  bolometric  magnitude  »?b0i  = 
3.08,  with  an  amplitude  Amm  =  0.49.  Over  the  last  100  yr,  the 
variability  period  of  S  Ori  has  been  seen  to  randomly  vary  be¬ 
tween  about  400  and  450  days  (Bedding  et  al.  1999;  Whitelock 
et  al.  2000;  Merchan  Benitez  &  Jurado  Vargas  2002).  For 
our  study,  we  adopt  a  stellar  period  of  420  days  and  a  time 
of  last  maximum  visual  brightness  of  JD  2,452,33 1 ,  as  given  by 
Pojmanski  (2002).  This  is  consistent  within  approximately 
±5  days  with  an  inspection  of  the  most  recent  AAVSO  and 
AFOEV  data.  Van  Belle  et  al.  (1996)  measured  the  near-infrared 
A'-band  uniform  disk  ( UD)  diameter  of  S  Ori  at  variability  phase 
0.56  to  be  10.54  ±  0.68  mas  and  derived  a  Rosseland  mean 
diameter  of  1 1 .70  ±  0.75  mas.  We  adopt  the  distance  to  S  Ori  to 
be  422  ±37  pc,  as  given  by  van  Belle  et  al.  (1996),  based  on 
measurements  by  Wyatt  &  Cahn  (1983)  and  Young  (1995).  This 
is  consistent  with  the  value  of  430  pc  from  Whitelock  et  al. 
(2000)  based  on  their  K  period-luminosity  relation.  S  Ori  ex¬ 
hibits  SiO  and  OH  masers  ( Benson  et  al.  1990),  but  no  detection 
of  any  H20  maser  emission  has  been  reported.  Sloan  &  Price 
(1998)  report  a  relatively  low  dust  emission  coefficient  (DEC), 
i.e.,  the  total  emission  of  the  dust  to  the  total  emission  of  the  star 
in  the  wavelength  range  from  7.7  to  14.0  /im,  of  0.24  for  S  Ori 
(for  comparison,  DEC  =  0.46  for  o  Cet). 

3.  OBSERVATIONS  AND  REDUCTION 

3.1.  Near-Infrared  K-Band  Interferometry 

We  obtained  near-infrared  A'-band  interferometric  measure¬ 
ments  of  S  Ori  using  the  ESO  Very  Large  Telescope  Interfer¬ 
ometer  (VLTI)  equipped  with  the  commissioning  instrument 
VINCI  and  the  two  VLTI  test  siderostats  (effective  apertures 
~40  cm).  VLTI  stations  D1  and  B3,  fomiing  an  unprojected 
baseline  length  of  24  m,  were  used  for  13  nights  between  2003 
January  25  (JD  2,452,665)  and  2003  March  3 1  (JD  2,452,730). 
The  S  Ori  visual  variability  phase  ranged  from  0.80  to  0.95 
for  these  dates.  Available  data  taken  in  2002  December  using 
an  8  m  ground  baseline  and  in  2003  September  using  a  16  m 
ground  baseline  were  not  used,  because  S  Ori  was  only  margin¬ 
ally  resolved  with  these  shorter  baselines  (|  V\  ^  0.9).  A  recent 
general  description  of  the  VLTI  can  be  found  in  Glindemann 
et  al.  (2003)  and  references  therein.  During  each  observation 
night,  several  series  of  typically  500  interferogram  scans  were 


TABLE  1 

Properties  of  the  Observed  VLTI-VINCI  Calibration 
Stars  from  Borde  et  al.  (2002) 


Star 

WUD 

(mas) 

±0) 

(mas) 

Teff 

(K) 

18  Mon . 

1.86 

0.023 

4656 

31  Ori . 

3.56 

0.057 

4046 

58  Hya . 

3.12 

0.035 

4318 

a  t  Ala . 

5.94 

0.016 

9900 

8  Lep . 

2.56 

0.041 

4656 

e  Lep . 

5.91 

0.064 

4046 

HD  112213 . 

3.15 

0.036 

3690 

HD  132833  . 

3.06 

0.034 

3690 

HR  2305 . 

1.76 

0.031 

4256 

HR  2311 . 

2.43 

0.040 

4046 

HR  3803 . 

6.93 

0.079 

4046 

L’  Sgr . 

3.28 

0.036 

4508 

i/2  (  Ala . 

2.38 

0.026 

4497 

T  Sgr . 

3.83 

0.043 

4444 

recorded  on  S  Ori  as  well  as  on  several  calibration  stars.  The 
calibration  stars  used  and  their  adopted  properties  are  listed  in 
Table  1.  For  S  Ori,  an  effective  temperature  of  2500  K  was 
used,  which  is  consistent  with  the  value  of2313±110K  ob¬ 
tained  by  van  Belle  et  al.  (1996)  for  phase  0.56.  We  note  that 
variations  of  the  effective  temperatures  up  to  ±500  K  result  in 
variations  of  the  squared  visibility  amplitude  of  less  than  0.3%. 
The  scan  length  for  the  S  Ori  and  all  calibration  stars’  obser¬ 
vations  was  230  /im,  the  scan  speed  was  616  /im  s  ,  and  the 
fringe  frequency,  corresponding  to  the  time  to  scan  one  inter¬ 
ferometric  fringe,  was  289.5  Hz. 

Mean  coherence  factors  were  obtained  for  each  series  of 
interferograms  using  the  VINCI  data  reduction  software,  ver¬ 
sion  3.0,  as  described  by  Kervella  et  al.  (2004),  employing  the 
results  based  on  wavelet  transforms.  Calibrated  squared  visi¬ 
bility  values  for  S  Ori  were  obtained  by  calibration  of  the  mean 
coherence  factors  as  described  in  Wittkowski  et  al.  (2004),  with 
a  time  kernel  of  5  hr  to  convolve  the  measured  transfer  func¬ 
tion.  The  errors  of  the  calibrated  S  Ori  squared  visibility  values 
include  the  scatter  of  the  single-scan  coherence  factors,  the 
adopted  errors  of  the  diameters  of  the  calibration  stars,  and  the 
observed  variation  of  the  transfer  function  during  each  night. 
The  calibrated  visibility  values  are  listed  in  Table  2  and  are  also 
available  in  electronic  form  from  the  authors  upon  request. 
Table  2  also  lists  the  date  and  time  of  observation,  the  spatial 
frequency,  the  position  angle  of  the  projected  baseline  (east  of 
north),  the  calibrated  squared  visibility  value  V2  and  its  error 
( 7Vi ,  as  well  as  the  number  of  processed  interferograms.  The 
effective  wavelength  is  2.19  /im.  All  observations  were  per¬ 
formed  at  similar  angles  of  the  projected  baseline  ranging  from 
72°  to  74°  east  of  north.  The  squared  visibility  values  for  S  Ori 
were  grouped  into  four  bins  of  five  nights  each  (bin  width  ~1% 
of  the  visual  variability  period  of  420  days). 

We  characterize  our  measurements  with  a  best-fitting  UD 
diameter  for  each  of  our  four  epochs.  Our  uniform  disk  model 
takes  into  account  the  broad  passband  (1.9-2. 5  /tm)  of  the 
VINCI  instrument,  as  described  in  Wittkowski  et  al.  (2004). 
Figure  1  shows  the  obtained  squared  visibility  amplitudes  of  S 
Ori  together  with  the  best-fitting  uniform  disk  models  for  each 
of  our  four  VLTI-VINCI  epochs.  The  larger  scatter  of  the  vis¬ 
ibility  data  for  the  March  16-21  epoch,  as  compared  to  the 
other  epochs,  can  be  explained  by  a  larger  time  difference  of 
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TABLE  2 

Calibrated  VLTI-VINCI  Visibility  Values 


Date 

Time 

Spatial  Frequency 
(arcsec-1) 

Baseline  Angle 
(deg  east  of  north) 

V2 

C 7y2 

Number  of  Processed  Interferograms 

2003  Jan  26 . 

01:56:48 

52.88 

72.23 

4.447E- 1 

1.053E— 2 

479 

2003  Jan  26 . 

02:02:42 

52.98 

72.38 

4.431E-1 

1.056E— 2 

458 

2003  Jan  26 . 

02:09:00 

53.05 

72.53 

4.358E-1 

1.062E— 2 

416 

2003  Jan  26 . 

02:47:04 

52.68 

73.19 

4.465E—  1 

1.093E— 2 

453 

2003  Jan  26 . 

02:53:21 

52.49 

73.26 

4.462E- 1 

1.101E— 2 

436 

2003  Jan  26 . 

02:59:39 

52.26 

73.31 

4.455E—  1 

1.106E— 2 

462 

2003  Jan  26 . 

03:11:06 

51.76 

73.38 

4.689E—  1 

1.194E— 2 

432 

2003  Jan  26 . 

03:17:07 

51.45 

73.40 

4.576E— 1 

1.208E— 2 

412 

2003  Jan  26 . 

03:23:36 

51.08 

73.41 

4.639E— 1 

1.495E— 2 

377 

2003  Jan  26 . 

04:11:15 

47.21 

73.03 

5.237E— 1 

2.989E— 2 

256 

2003  Jan  26 . 

04:17:31 

46.56 

72.92 

5.417E—  1 

3.260E— 2 

250 

2003  Jan  26 . 

04:36:44 

44.37 

72.45 

5.615E—  1 

2.009E— 2 

241 

2003  Jan  26 . 

04:43:26 

43.55 

72.25 

5.695E-1 

1.948E— 2 

210 

2003  Jan  31 . 

02:40:32 

52.24 

73.32 

4.838E— 1 

4.282E— 2 

145 

2003  Jan  31 . 

02:51:39 

51.75 

73.38 

4.662E—  1 

4.361E— 3 

479 

2003  Jan  31 . 

02:57:52 

51.43 

73.40 

4.753E— 1 

4.382E— 3 

485 

2003  Jan  31 . 

03:49:26 

47.43 

73.07 

5.284E—  1 

7.397E— 3 

365 

2003  Feb  5 . 

03:53:23 

44.86 

72.56 

5.598E—  1 

1.807E— 2 

364 

2003  Feb  5 . 

04:05:16 

43.40 

72.21 

5.980E-1 

2.105E— 2 

232 

2003  Feb  6 . 

03:26:20 

47.38 

73.06 

5.323E— 1 

1.217E— 2 

471 

2003  Feb  6 . 

03:58:48 

43.72 

72.29 

5.949E—  1 

1.483E— 2 

388 

2003  Feb  10 . 

01:08:45 

53.04 

72.50 

5.110E— 1 

4.295E— 2 

159 

2003  Feb  10 . 

01:28:27 

53.03 

72.91 

4.723E— 1 

2.416E— 2 

427 

2003  Feb  10 . 

02:17:13 

51.50 

73.40 

4.962E—  1 

2.543E— 2 

454 

2003  Feb  10 . 

02:23:25 

51.15 

73.41 

4.950E-1 

2.536E— 2 

456 

2003  Feb  10 . 

02:29:57 

50.74 

73.41 

5.159E— 1 

2.668E— 2 

443 

2003  Feb  10 . 

03:24:02 

45.96 

72.80 

5.803E—  1 

4.593E— 2 

250 

2003  Mar  17 . 

00:26:00 

49.77 

73.36 

5.115E— 1 

7.606E— 3 

323 

2003  Mar  17 . 

00:33:21 

49.18 

73.30 

5. 12  IE —  1 

1.296E— 2 

245 

2003  Mar  17 . 

00:42:06 

48.42 

73.21 

5.156E— 1 

1.454E— 2 

192 

2003  Mar  18 . 

00:12:28 

50.47 

73.40 

5.041E—  1 

7.892E— 3 

459 

2003  Mar  18 . 

00:25:52 

49.47 

73.33 

5.173E— 1 

9.933E— 3 

309 

2003  Mar  18 . 

00:33:16 

48.85 

73.27 

5.271E—  1 

1.035E— 2 

345 

2003  Mar  21 . 

00:43:22 

46.75 

72.95 

5.772E— 1 

9.288E— 3 

385 

2003  Mar  21 . 

00:48:38 

46.19 

72.85 

6.055E—  1 

2.220E— 2 

42 

2003  Mar  21 . 

01:06:55 

44.07 

72.38 

6.018E— 1 

1.650E— 2 

323 

2003  Mar  22 . 

00:34:15 

47.29 

73.04 

5.797E—  1 

4.487E— 3 

454 

2003  Mar  22 . 

00:41:19 

46.55 

72.92 

5.837E— 1 

4.348E— 3 

471 

2003  Mar  22 . 

00:48:31 

45.77 

72.76 

5.960E—  1 

4.632E— 3 

461 

2003  Mar  27 . 

23:51:52 

49.03 

73.29 

5.253E—  1 

1.073E— 2 

460 

2003  Mar  27 . 

23:59:06 

48.39 

73.21 

5.368E— 1 

1.139E— 2 

389 

2003  Mar  28 . 

00:06:22 

47.71 

73.11 

5.500E—  1 

1.248E— 2 

322 

2003  Mar  28 . 

23:43:31 

49.40 

73.32 

5.304E—  1 

1.179E— 2 

400 

2003  Mar  28 . 

23:50:53 

48.78 

73.26 

5.279E— 1 

1.467E— 2 

348 

2003  Mar  28 . 

23:58:29 

48.09 

73.17 

5.558E-1 

2.190E— 2 

272 

2003  Mar  29 . 

23:34:06 

49.83 

73.36 

5.201E—  1 

9.993E— 3 

447 

2003  Mar  29 . 

23:41:37 

49.23 

73.31 

5.238E— 1 

1.175E— 2 

352 

2003  Mar  29 . 

23:49:13 

48.58 

73.23 

5.424E—  1 

1.325E— 2 

335 

2003  Mar  30 . 

00:31:28 

44.08 

72.38 

6.009E- 1 

1.13  IE — 2 

459 

2003  Mar  3 1 . 

23:37:44 

48.89 

73.27 

5.267E—  1 

2.318E— 2 

342 

2003  Apr  1 . 

00:03:01 

46.45 

72.89 

5.604E—  1 

1.989E— 2 

402 

the  calibration  stars.  Table  3  lists  the  observation  dates  and  the 
determined  A-band  UD  diameters  and  their  associated  formal 
errors.  Additional  calibration  uncertainties  are  estimated  to  be 
~0. 1  mas,  based  on  an  analysis  of  subvolumes  of  the  data  and 
differing  methods  of  interpolating  the  transfer  function.  The 
addition  of  a  possible  circumstellar  dust  shell  is  not  expected  to 
affect  the  A'-band  UD  diameter  significantly  (see,  for  instance, 
discussions  in  Woodruff  et  al.  2004;  Ohnaka  et  al.  2005).  In 
particular,  Ohnaka  et  al.  (2005)  modeled  an  optically  thin  dust 


disk  for  a  similar  Mira  star,  RR  Sco  (DEC  =  0.21,  compared  to 
the  value  of  0.24  for  S  Ori),  for  the  purpose  of  comparison  to 
mid-infrared  interferometric  data.  At  a  wavelength  of  2.2  /im 
the  scattered  and  thermal  emission  from  this  dust  disk  is  lower 
by  a  factor  of  ~1000  than  the  attenuated  star  light,  and  the 
effect  on  the  squared  visibility  values  at  V 2  ~  0.5  is  less  than 
0.2%  (K.  Ohnaka  2004,  private  communication).  This  would 
lead  to  an  effect  on  the  UD  diameters  in  Table  3  of  less  than 
0.02  mas.  Discussed  in  §  4.1  below  is  the  relationship  between 
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Fig.  1. — Near-infrared  A^-band  squared  visibility  amplitudes  of  S  Ori  and  best-fitting  uniform  disk  models  for  our  four  VLTI- VINCI  epochs. 


a  K- band  UD  diameter  and  a  physically  more  meaningful  an¬ 
gular  size  of  the  stellar  photosphere. 

3.2.  SiO  Maser  Observations 

We  simultaneously  observed  the  v  =  1,  J  =  1-0  (43.1  GHz) 
and  v  =  2,  J  =  1-0  (42.8  GHz)  SiO  maser  transitions  toward 
S  Ori  (a  =  05h29m00?9,  6  =  -04°41'3277,  J2000.0)  on  2002 
December  29  for  6  hr  starting  at  02:43  UT  (JD  2,452,636.6, 
visual  variability  phase  0.73).  S  Ori  and  a  continuum  calibrator 
(0359+509)  were  observed  using  the  10  stations  of  the  VLBA. 
The  VLBA  is  operated  by  the  National  Radio  Astronomy  Ob¬ 
servatory  (NRAO).1  Reference  frequencies  of  43. 122080  and 
42.820587  GHz  were  used  for  the  v  =  1  and  v  =  2  SiO  transi¬ 
tions,  respectively.  Data  were  recorded  in  dual  circular  polar¬ 
ization  using  two  8  MHz  (56.1  km  s  1 )  bands  centered  on  the 
local  standard  of  rest  (LSR)  velocity  of  18.0  km  s  1 .  System 
temperatures  and  point-source  sensitivities  were  on  the  order 
of  ~  1 50  K  and  ^11  Jy  K-1,  respectively. 

1  The  National  Radio  Astronomy  Observatory  is  a  facility  of  the  National 
Science  Foundation  operated  under  cooperative  agreement  by  Associated 
Universities,  Inc. 


The  data  were  correlated  at  the  VLBA  correlator  operated  by 
NRAO  in  Socorro,  New  Mexico.  Auto-  and  cross-correlation 
spectra  consisting  of  256  channels  with  channel  spacings  of 
31.25  kHz  (~0.2  km  s  1 )  were  produced  by  the  correlator. 
Calibration  was  performed  using  the  Astronomical  Image  Pro¬ 
cessing  System  (AIPS)  maintained  by  NRAO.  The  total  in¬ 
tensity  data  were  calibrated  in  accordance  with  the  procedures 
outlined  in  Diamond  (1989).  The  bandpass  response  was  de¬ 
termined  from  scans  on  the  continuum  calibrator  and  used  to 
correct  the  target  source  data.  The  time-dependent  gains  of  all 
antennas  relative  to  a  reference  antenna  were  determined  by 
fitting  a  total  power  spectrum  (from  the  reference  antenna  with 
the  target  source  at  a  high  elevation)  to  the  total  power  spectrum 
of  each  antenna.  The  absolute  flux  density  scale  was  established 
by  scaling  these  gains  by  the  system  temperature  and  gain  of  the 
reference  antenna.  Errors  in  the  gain  and  pointing  of  the  refer¬ 
ence  antenna  and  the  atmospheric  opacity  contribute  to  the  error 
in  the  absolute  amplitude  calibration,  which  is  accurate  to  about 
15%-20%. 

To  correct  any  instrumental  delay,  a  fringe  fit  was  performed 
on  the  continuum  calibrator  scans,  and  residual  group  delays  for 
each  antenna  were  determined.  Variations  in  the  residual  delays 
ranged  from  2  to  4  ns  resulting  in  phase  errors  of  no  more  than 


TABLE  3 

VLTI- VINCI  Observations  and  Derived  A>Band  UD  Diameters  with  Associated  Formal  Errors 


Gud  0(6) 

Date  Mean  JD  Number  of  Series  Mean  Phase  (mas)  (mas)  xl 


Jan  25-30 .  2,452,666.8  17  0.80  10.52  0.03  0.50 

Feb  04-09 .  2,452,678.8  10  0.83  10.33  0.07  0.55 

Mar  16-21 .  2,452,718.0  12  0.92  9.98  0.06  3.91 

Mar  27-31 .  2,452,728.1  12  0.95  10.16  0.03  0.13 


Note. — The  observation  dates  were  grouped  into  four  bins  of  5  days  each. 
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Fig.  2. — Total  intensity  image  of  the  v  =  1,  /  =  1—0  (red)  and  v  =  2,  J  = 
1—0  (blue)  SiO  maser  emission  toward  S  Ori.  The  images  represent  maximum 
pixel  values  over  the  LSR  velocity  range  from  +10.1  to  +25.4  km  s'  plotted 
as  contours.  Contour  levels  are  1,  2,  4,  8,  16,  32,  and  64  times  the  3  a  off- 
source  noise  of  0.1  mjy  beam-1.  Peak  flux  density  is  4.36  Jy  beam-1  for  the 
v  =  1,  J  =  1—0  line  and  5.35  Jy  beam-1  for  the  v  —  2,  ./  —  1—0  transition. 
Synthesized  beam  sizes  are  0.58  x  0.17  mas  at  a  position  angle  of  —  1 5.4  for 
v  =  1 ,  J  =  1—0  and  0.52  x  0. 1 7  mas  at  a  position  angle  of  —  1 5?7  for  v  =  2, 
J  =  1-0. 

1?5— 3°  across  the  8  MHz  band.  Residual  fringe  rates  were  ob¬ 
tained  by  fringe-fitting  a  strong  reference  feature  in  the  spec¬ 
trum  of  each  maser  transition.  For  both  transitions  we  used  the 
same  channel  at  a  velocity  Flsr=  16.2  km  s~*.  The  resulting 
fringe  rate  solutions  were  applied  to  all  channels  in  each  spec¬ 
trum,  respectively.  An  iterative  self-calibration  and  imaging 
procedure  was  then  performed  to  map  this  reference  channel 
for  each  transition.  The  resulting  residual  phase  and  amplitude 
corrections  from  the  reference  channels  at  42.8  and  43.1  GHz 
were  applied  to  all  channels  in  the  respective  bands. 

In  order  to  accurately  compare  the  distributions  of  the  two 
maser  transitions,  it  is  desirable  to  determine  a  common  spatial 
reference  point.  However,  after  the  fringe-fitting  step  to  de¬ 
termine  residual  fringe  rates,  all  absolute  position  information 
is  lost  for  the  VLBA  data.  We  accomplished  the  registration  of 
the  two  transitions  by  applying  the  fringe-fitted  solutions  from 
the  v  =  2,  42.8  GHz  transition  to  the  v  =  1,  43.1  GHz  data,  re¬ 
mapping  the  v  =  1  transition,  and  comparing  the  resulting  im¬ 
ages  with  those  produced  from  the  application  of  the  v  =  1 
calibration  itself.  Although  the  images  resulting  from  the  ap¬ 
plication  of  the  42.8  GHz  fringe  fit  to  the  43.1  GHz  data  were 
of  poorer  quality,  we  were  still  able  to  determine  an  offset  be¬ 
tween  the  two  sets  of  images.  The  offsets  in  right  ascension  and 
declination  computed  from  images  of  two  different  spectral 
channels  were  the  same  to  within  0.03  mas.  Subsequent  phase 
self-calibration  and  imaging  were  performed  for  the  v  =  1  SiO 
data  using  these  computed  offsets,  resulting  in  spatially  aligned 
v  =  1  and  v  =  2  image  cubes.  This  same  procedure  was  ap¬ 
plied  to  the  continuum  calibrator  source,  0359+509,  and  the 
positions  derived  agreed  to  within  0.4  mas,  thus  providing 
an  estimate  of  the  error  in  the  registration  of  the  two  maser 
distributions. 


Final  images  of  the  SiO  maser  emission  consisting  of  1024  x 
1024  pixels  (~51  x  51  mas)  were  generated  using  synthesized 
beams  of  0.58x0. 17  mas  and  0.52  x  0.17  mas  for  the  v  =  1  and 
v  =  2  transitions,  respectively.  Images  were  produced  for  spec¬ 
tral  channels  from  4.7  to  30.7  km  s  ,  forming  image  cubes  of 
120  planes.  Off-source  rms  noise  in  the  images  ranged  from  7  to 
18  mjy.  Figure  2  shows  the  total  intensity  contour  maps  of  the 
v  =  1,  43.1  GHz  {red)  and  v  =  2,  42.8  GHz  {blue)  SiO  maser 
emission  toward  S  Ori.  The  contours  represent  the  maximum 
pixel  value  in  each  image  cube  over  the  LSR  velocity  range 
from+10.1  to  +25.4kms~1.  Subsequent  analysis  of  the  image 
data  is  described  in  detail  in  §  4.2. 

4.  RESULTS  AND  DISCUSSION 
4.1.  The  Photospheric  Diameter 

Figure  3  and  Table  3  show  our  near-infrared  A-band  UD 
diameters  of  S  Ori  determined  from  our  four  VLTI- VINCI 
epochs  as  a  function  of  observation  date  and  visual  variability 
phase,  as  discussed  in  §  3.1.  The  A'-band  UD  diameter  de¬ 
creases  almost  linearly  from  ~10.5  mas  at  visual  variability 
phase  0.80  to  ~10.2  mas  at  phase  0.95,  i.e.,  by  ^3%.  It  is  not 
clear  whether  the  local  diameter  minimum  of  ~ 1 0.0  mas  at 
phase  0.92  in  Table  3  and  Figure  3  is  caused  by  an  additional 
systematic  calibration  error,  as  discussed  in  §  3.1,  or  by  a  real 
minimum  of  the  stellar  size.  The  minimum  of  a  Mira  A’-band 
UD  diameter  is  predicted  to  lie  at  visual  variability  phase  ~0.9 
according  to  model  predictions  by  Ireland  et  al.  (2004).  Our  co¬ 
ordinated  VLBA  measurement  at  visual  variability  phase  0.73 
occurred  about  one  month  before  the  first  VLTI  epoch  at  phase 
0.8  and  is  hence  close  to,  but  not  exactly  contemporaneous  to, 
the  VLTI  epochs.  Both  observations  of  the  diameter  change  of 
the  Mira  star  S  Lac  by  Thompson  et  al.  (2002b)  and  Mira  star 
model  predictions  by  Ireland  et  al.  (2004)  show  an  almost 
linear  change  of  the  A'-band  UD  diameter  between  variability 
phases  0.7  and  0.9,  i.e.,  in  the  preminimum  part  of  the  diameter 
curve  and  premaximum  part  of  the  visual  light  curve.  Thus,  we 
can  linearly  extrapolate  our  A-band  UD  diameter  data  at  phases 
0.8-0.95  to  the  VLBA  phase  0.73.  For  phase  0.73,  we  deter¬ 
mine  an  extrapolated  A-band  UD  diameter  of  0yD(VLBA 
epoch,  phase  =  0.73)  =  10.8  ±  0.3  mas.  This  value  and  its 


JD -  2450000 

2650  2700  2750  2800  2850  2900  2950 


Fig.  3. — Obtained  S  Ori  K- band  UD  diameter  as  a  function  of  date  and 
variability  phase.  The  dashed  lines  denote  the  best-fitting  linear  functions  to 
all  of  our  values  and  to  the  two  closest  points  to  our  VLBA  epoch.  Our  VLBA 
and  VLTI  variability  phases  are  indicated  by  arrows. 
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Fig.  4. — LOS  velocity  structure  of  the  v  —  1,  J  —  1—0  SiO  maser  emission 
toward  S  Ori.  The  top  panel  shows  the  spectrum  formed  by  plotting  maser 
intensity  vs.  velocity,  color  coded  in  2  km  s_1  velocity  increments  from  red- 
ward  (left)  to  blueward  (right).  The  solid  line  in  the  top  panel  represents  the 
scalar-averaged  cross-power  spectrum  averaged  over  all  of  the  VLBA  anten¬ 
nas.  The  bottom  panel  plots  the  spatial  and  velocity  distribution  of  the  masers. 
The  color  of  each  point  represents  the  corresponding  velocity  bin  in  the 
spectrum,  and  the  size  of  each  point  is  proportional  to  the  logarithm  of  the  flux 
density.  Errors  in  the  positions  of  the  features  are  smaller  than  the  data  points. 
The  dashed  circle  is  based  on  the  mean  angular  distance  of  the  SiO  masers  from 
the  center  of  the  distribution.  The  colored  circle  in  the  center  shows  the  angular 
size  of  the  photosphere  as  determined  from  our  VLTI  K- band  measurements. 


Fig.  5. — LOS  velocity  structure  of  the  v  —  2,  J  —  1—0  SiO  maser  emission 
toward  S  Ori.  The  top  panel  shows  the  spectrum  formed  by  plotting  maser 
intensity  versus  velocity,  color  coded  in  2  km  s-1  velocity  increments  from 
redward  (left)  to  blueward  (right).  The  solid  line  in  the  top  panel  represents  the 
scalar-averaged  cross-power  spectrum  averaged  over  all  of  the  VLBA  anten¬ 
nas.  The  bottom  panel  plots  the  spatial  and  velocity  distribution  of  the  masers. 
The  color  of  each  point  represents  the  corresponding  velocity  bin  in  the 
spectrum,  and  the  size  of  each  point  is  proportional  to  the  logarithm  of  the  flux 
density.  Errors  in  the  positions  of  the  features  are  smaller  than  the  data  points. 
The  dashed  circle  is  based  on  the  mean  angular  distance  of  the  SiO  masers  from 
the  center  of  the  distribution.  The  colored  circle  in  the  center  shows  the  angular 
size  of  the  photosphere  as  determined  from  our  VLTI  A^-band  measurements. 


error  correspond  to  the  mean  and  difference  of  the  two  linear 
extrapolations  shown  in  Figure  3,  which  used  all  four  points 
and  only  the  two  closest  points  to  the  VLBA  epoch,  respec¬ 
tively.  This  UD  diameter  corresponds  to  a  linear  radius  R^D  = 
2.3  ±  0.5  AU  or  =  490  ±115  R.,  with  our  assumed  dis¬ 
tance  to  S  Ori. 

A  UD  intensity  profile  is  often  not  an  ideal  representation  of 
the  true  near-infrared  center-to-limb  intensity  variation  (CLV) 
of  cool  giants  in  general  and  of  Mira  stars  in  particular.  The 
CLV  of  cool  giants  has  been  studied  by,  e.g.,  Quirrenbach  et  al. 
(1996),  Bums  et  al.  (1997),  Scholz  (1998),  Flofmann  et  al. 
(1998),  Wittkowski  et  al.  (2001,  2004),  Ireland  et  al.  (2004), 
and  Woodruff  et  al.  (2004).  UD  diameters  of  nonvariable  giants 
obtained  from  visibility  measurements  in  the  first  lobe  are  usu¬ 
ally  transformed  into  physically  more  meaningful  Rosseland 
angular  diameters  using  correction  factors  determined  from 
atmosphere  models  in  the  literature  (e.g.,  Claret  2003).  For 
the  case  of  Mira  variable  stars,  few  atmosphere  models  are 
available.  We  used  the  hydrodynamic  atmosphere  models  from 
Hofmann  et  al.  (1998),  Tej  et  al.  (2003),  and  Ireland  et  al. 
(2004)  to  estimate  the  relationship  between  the  W-band  UD  di¬ 
ameter  and  the  continuum  diameter.  While  the  broadband  UD 


diameter  is  affected  by  various  molecular  bands,  the  contin¬ 
uum  diameter  is  a  better  estimate  of  the  real  photospheric 
stellar  size.  These  models  have  been  constructed  for  the  pro¬ 
totype  Mira  stars  o  Cet  and  R  Leo  and  have  been  compared  to 
observations  of  several  Mira  stars  by,  for  instance,  Hofmann 
et  al.  (2001,  2002)  and  Woodruff  et  al.  (2004).  The  general 
model  results  are  not  expected  to  be  dramatically  different 
for  other  Mira  stars  such  as  S  Ori  (M.  Scholz  2004,  private 
communication).  The  model  predicted  difference  between  the 
continuum  and  V-band  UD  diameters  is  relatively  low  in  the 
premaximum  region  of  the  visual  variability  curve,  as  in  the  case 
of  our  observations.  At  this  phase  of  0.73,  the  continuum  di¬ 
ameter  may  be  smaller  than  the  W-band  UD  diameter  by  about 
15%  (Ireland  et  al.  2004).  With  this  assumption,  the  continuum 
photospheric  diameter  for  the  epoch  of  our  VLBA  observation 
would  be  0phot(VLBA  epoch,  phase  =  0.73)  «  9.2  mas.  This 
angular  radius  corresponds  to  a  photospheric  radius  of  Rphot  ~ 
420  Rq  or  /x’phot «  1.9  AU  with  our  assumed  distance  to  S  Ori. 
With  a  bolometric  flux  of  m\,0\  =  3.1  ±  0.3  (Whitelock  et  al. 
2000),  our  value  for  0phOt  corresponds  to  an  effective  tempera¬ 
ture  TeS  ~  2670  K. 
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4.2.  The  v  =  1  and  v  =  2  SiO  Maser  Emission 

Our  total  intensity  images  of  the  43.1  and  42.8  GHz  SiO 
maser  emission  (Fig.  2)  show  a  typical  clumpy  distribution 
of  the  maser  spots  within  a  ringlike  structure.  For  the  v  =  2, 
42.8  GHz  transition  the  ring  is  relatively  sparse  with  nearly  all 
of  the  masers  concentrated  to  the  northwest  side  of  the  shell 
and  a  few  features  to  the  southeast.  The  v=  1,  43.1  GHz  ma¬ 
sers,  however,  form  a  more  typical  ringlike  structure  often  seen 
for  other  stars  with  SiO  masers.  Like  the  v  =  2  masers,  the 
v  =  1  SiO  also  has  a  higher  concentration  of  features  on  the 
northwest  side  of  the  shell.  The  v  =  1  maser  ring  appears  to  be 
symmetrical,  forming  a  nearly  circular  shell.  The  distribution 
of  the  v  =  2  SiO  masers  is  too  sparse  to  allow  us  to  comment 
on  any  symmetry  for  this  transition. 

In  order  to  identify  and  parameterize  maser  components, 
two-dimensional  Gaussian  functions  were  fitted  to  the  emission 
in  each  spectral  (velocity)  channel  using  the  A1PS  task  SAD. 
Image  quality  was  assessed  using  the  off-source  rms  noise  and 
the  deepest  negative  pixel  in  the  image.  A  cutoff  flux  density 
was  conservatively  set  to  the  greater  of  8  errms  or  the  absolute 
value  of  the  deepest  negative  pixel  in  the  plane.  Features  with 
flux  densities  greater  than  this  cutoff  were  fit  with  Gaussians 
to  determine  component  parameters.  Errors  in  the  positions  in 
right  ascension  and  declination  of  identified  features  were  com¬ 
puted  using  the  fitted  source  size  divided  by  twice  the  signal- 
to-noise  ratio  (S/N)  in  the  image  and  ranged  from  1  /ras  for 
features  with  high  S/N  to  50  /ras  for  features  with  lower  S/N. 

Since  the  ^0.2  km  s~ 1  channel  spacing  is  sufficient  to  resolve 
the  masers  spectrally,  features  typically  appear  in  multiple  ad¬ 
jacent  spectral  channels.  Positions  in  right  ascension  and  dec¬ 
lination  and  center  velocities  for  the  masers  were  determined 
using  a  flux-density-squared  weighted  average  for  features 
identified  in  two  or  more  adjacent  channels  with  a  spatial  co¬ 
incidence  of  0.2  mas  (~2/3  of  the  geometric  mean  of  the  syn¬ 
thesized  beam).  The  flux  assigned  to  the  maser  averages  was  the 
maximum  single-channel  flux  density.  The  maser  components 
identified  using  this  procedure  are  represented  by  the  circles  in 
Figures  4  and  5.  In  the  figures,  point  sizes  are  proportional  to 
the  logarithm  of  the  fitted  flux  density,  which  ranged  from  0.2  to 
4.3  Jy  for  Figure  4  and  0.1  to  5.3  Jy  for  Figure  5.  The  line-of- 
sight  (LOS)  velocity  information  for  the  masers  is  also  repre¬ 
sented  in  Figures  4  and  5 .  The  top  panel  of  each  figure  shows  the 
spectrum  of  SiO  maser  emission  ranging  from  10  to  25  km  s_1 
color  coded  by  LOS  velocity  in  increments  of  2  km  s_1.  The 
bottom  panels  of  Figures  4  and  5  show  the  spatial  distribution 
of  the  SiO  masers  plotted  with  the  same  velocity  color  coding  as 
in  the  top  panel  and  with  the  color  of  the  maser  representing  its 
corresponding  velocity  range  in  the  spectrum.  Comparing  the 
maser  component  maps  with  the  total  intensity  contour  maps, 
we  see  that  the  identified  features  accurately  represent  the  emis¬ 
sion  summed  over  all  velocity  channels  in  the  image  cube.  From 
the  component  maps,  there  does  not  appear  to  be  any  coherent 
velocity  structure  indicative  of  global  expansion /infall  or  ro¬ 
tation.  We  do  note  that  a  group  of  components  on  the  western 
side  of  the  shell  shows  a  velocity  gradient  with  velocity  de¬ 
creasing  with  increasing  distance  from  the  star.  Such  velocity 
gradients  in  the  SiO  maser  emission  have  been  observed  pre¬ 
viously  (e.g.,  Boboltz  &  Marvel  2000;  Hollis  et  al.  2001). 

To  compare  the  size  of  the  SiO  maser  distribution  to  the 
photospheric  diameter  of  the  star  as  measured  by  the  near- 
infrared  K- band  interferometry,  we  determined  the  average 
distance  of  the  masers  from  the  center  of  the  distribution.  To  ac¬ 
complish  this,  we  first  determined  the  center  of  the  distribution 


by  performing  a  least-squares  fit  of  a  circle  to  the  combined 
v  =  1  and  v  =  2  maser  component  data.  This  fit  produced  a 
common  center  from  which  we  computed  the  mean  maser  an¬ 
gular  distance  Fsio  and  the  standard  deviation  for  each  transi¬ 
tion  independently.  The  mean  angular  distances  from  the  center 
for  the  observed  SiO  masers  at  v  =  1 ,  43.1  GHz  and  v  =  2, 
42.8  GHz  are  Fs;o  =  9.4  and  8.8  mas,  respectively.  These  dis¬ 
tances  are  indicated  by  a  dashed  circle  and  are  listed  in  the 
bottom  panels  of  Figures  4  and  5.  Standard  deviations  of  the 
distances  are  <7sio  =1.4  and  1.7  mas  for  the  v  =  1  and  v  =  2 
transitions,  respectively.  The  standard  deviations  provide  an 
indication  of  the  thickness  of  the  shell  and  are  likely  dominated 
by  the  features  on  the  western  side  of  the  shell  with  a  wide  range 
of  distances  from  the  center.  The  mean  angular  distances  that 
we  derive  are  consistent  with  least-squares  circle  fits  to  each 
distribution  independently  using  the  previously  mentioned  com¬ 
mon  center.  The  least-squares  fits  to  each  transition  have  the 
added  assumption  that  the  distribution  is  circular;  thus,  we  have 
chosen  to  report  only  the  mean  distance  from  center.  The  com¬ 
puted  angular  distances  translate  to  linear  shell  distances  from 
the  center  of  4.0  ±  0.6  AU  and  3.7  ±  0.6  AU  at  the  assumed 
distance  of  422  ±  37  pc.  The  standard  deviations  are  likewise 
0.6  ±0.1  AU  and  0.7  ±0.1  AU  for  the  v  =  1  and  v  =  2  masers, 
respectively. 

In  principle,  a  comparison  of  the  relative  spatial  locations  of 
the  v  =  1  and  v  =  2  masers  should  allow  us  to  comment  on  pos¬ 
sible  pumping  mechanisms  (i.e.,  collisional,  radiative,  or  com¬ 
bination),  as  in  Desmurs  et  al.  (2000).  However,  because  the 
v  =  2,  42.8  GHz  masers  are  mostly  confined  to  a  small  region 
on  the  northwest  side  of  the  shell,  it  is  difficult  to  unambigu¬ 
ously  determine  the  relative  shift  between  the  two  maser  shells. 
Although  the  mean  distance  from  center  that  we  determine  for 
the  v  =  2  masers  is  slightly  smaller  than  the  v  =  1  mean  dis¬ 
tance,  the  two  values  are  consistent  to  within  the  errors. 

4.3.  Comparison  o/SiO  Maser  Distribution 
and  Photospheric  Stellar  Diameter 

Figures  4  and  5  also  provide  a  comparison  of  the  distribution 
of  the  SiO  maser  spots  with  our  obtained  angular  size  of  the 
photospheric  disk  at  the  same  epoch.  The  average  distance  of 
the  maser  spots  from  the  center  of  their  distribution  for  the  43. 1 
and  42.8  GHz  transitions  (F§ ;o=  9.4  and  8.8  mas)  appears  at 
1.7 R*  and  1.6 R*,  respectively,  when  compared  to  our  estimate 
of  the  K- band  UD  diameter.  When  compared  to  our  estimate 
of  the  continuum  diameter,  these  values  are  slightly  larger  at 
2. OR*  and  1.9R*,  respectively.  The  widths  of  the  distributions 
are  approximately  0.3R*  for  the  masers  in  both  transitions.  As  a 
reference,  the  angular  photospheric  radius  is  listed  in  the  bot¬ 
tom  panels  of  Figures  4  and  5  and  is  indicated  by  a  colored 
circle  at  the  center  of  the  maser  distribution.  The  position  of  the 
near-infrared  stellar  disk  indicated  on  the  plots  is  assumed  to 
coincide  with  the  center  of  the  SiO  maser  distribution  dis¬ 
cussed  in  §  4.2.  The  true  location  of  the  star  relative  to  the 
masers  is  still  unknown.  We  note  that  the  above  estimates  are 
based  on  angular  sizes  derived  for  the  SiO  maser  emission  and 
the  near-infrared  photosphere  and  hence  are  independent  of  the 
distance  to  S  Ori.  This  result  is  consistent  with  theoretical 
estimates  by  Humphreys  et  al.  (2002)  based  on  a  stellar  hydro- 
dynamic  pulsation  model  combined  with  an  SiO  maser  model. 
They  obtain  values  between  1 .7R*  and  2. 1R*,  depending  on  the 
variability  phase. 

Our  result  is  consistent  with  an  SiO  maser  ring  radius  of 
approximately  2 R*,  as  determined  by  Cotton  et  al.  (2004)  for 


960 


BOBOLTZ  &  WITTKOWSKI 


Vol.  618 


several  Mira  variables.  Recently,  however,  Monnier  et  al.  (2004) 
updated  the  stellar  diameter  estimate  of  Greenhill  et  al.  ( 1 995)  for 
the  supergiant  VX  Sgr  and  found  that  the  SiO  masers  lie  at  a 
greater  distance  from  the  photosphere,  3.9R*,  rather  than  the 
1.3R*  determined  by  Greenhill  et  al.  (1995).  It  is  striking  that  the 
stellar  diameter  for  the  other  supergiant  observed  by  Monnier 
et  al.  (2004),  NML  Cyg,  also  indicates  an  SiO  maser  ring  dis¬ 
tance  of  ~4 R*  when  compared  to  the  SiO  maser  shell  diameter 
obtained  by  Boboltz  &  Marvel  (2000).  It  is  unclear  whether  these 
larger  distances  are  an  artifact  caused  by  the  noncontempora- 
neous  measurements  or  whether  they  indicate  an  inherent  differ¬ 
ence  between  Mira  variables  and  supergiants.  In  earlier  studies, 
measured  sizes  for  SiO  maser  shells  have  been  compared  to 
stellar  sizes  based  on  diameters  found  in  the  literature  (e.g., 
Boboltz  et  al.  1997;  Boboltz  &  Marvel  2000;  Hollis  et  al.  2001). 
These  diameters  are  often  widely  spaced  in  time  and  stellar  phase 
from  the  measurements  of  the  SiO  maser  shell.  Sometimes  the 
stellar  diameters  have  not  been  measured  at  all,  and  comparisons 
are  made  with  even  less  precise  diameter  estimates  based  on  the 
luminosity  of  the  star  (e.g.,  Diamond  et  al.  1994;  Colomer  et  al. 
1996;  Sanchez  Contreras  et  al.  2002).  These  examples  demon¬ 
strate  that  it  is  desirable  to  compare  SiO  maser  ring  diameters 
with  stellar  diameter  estimates  and  highlight  the  need  for  con¬ 
temporaneous  multiwavelength  observations. 

5.  SUMMARY 

Using  the  VLBA  and  VLT1-VINCI,  we  have  undertaken  a 
coordinated  multiwavelength  study  of  the  Mira  variable  S  Ori. 
The  VLBA  observations  resulted  in  the  first  interferometric 
images  of  the  v  =  1,  J  =  1-0,  43.1  GHz  and  v  =  2,  J  =  1-0, 
42.8  GHz  SiO  transitions.  These  images  show  that  the  masers 
lie  in  a  clumpy  ringlike  distribution.  Analysis  of  the  VLBA 
images  provided  mean  distances  from  the  center  of  9.4  mas 
(4.0  AU)  and  8.8  mas  (3.7  AU)  and  standard  deviations  of 
1.4  mas  (0.6  AU)  and  1.7  mas  (0.7  AU)  for  the  v  =  1  and  v  = 
2  SiO  masers,  respectively,  at  a  stellar  phase  of  0.73.  From 
near-infrared  A'- band  visibility  measurements  made  with  VLT1- 
VINCI,  we  derived  UD  diameters  that  decreased  over  time 
from  ~10.5  mas  at  stellar  phase  0.80  to  ~10.2  mas  at  phase 
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0.95.  After  extrapolating  the  UD  diameter  to  the  phase  of  our 
VLBA  observations  and  considering  a  correction  from  UD  di¬ 
ameter  to  continuum  photospheric  size,  we  determine  a  photo- 
spheric  diameter  of  0phot  ~  9.2  mas,  corresponding  to  a  linear 
radius  of  Rph0t  ~  1  -9  AU. 

Because  these  observations  of  S  Ori  were  closely  spaced  in 
time  in  a  coordinated  effort  between  the  two  instruments,  we 
are  able  to  relate  the  stellar  diameter  with  the  size  of  the  SiO 
maser  shell  without  the  uncertainty  caused  by  the  inherent  var¬ 
iability  of  the  star.  Thus,  we  can  conclusively  say  that  the  SiO 
masers  lie  relatively  close  to  the  stellar  photosphere  at  a  dis¬ 
tance  of  R*  at  the  time  of  our  observations. 

The  observations  presented  here  represent  the  first  in  a  series 
of  coordinated  VLBA/VLTI  experiments  to  study  long-period 
variable  stars.  The  present  VLTI  observations  were  obtained  dur¬ 
ing  the  commissioning  period  of  the  VLTI  with  the  commis¬ 
sioning  instrument  VINCI  and  the  small  (40  cm)  test  siderostats. 
Future  observations  will  use  the  scientific  VLTI  instruments 
in  the  mid-infrared  (MIDI)  and  near-infrared  (AMBER)  with 
the  8  m  unit  telescopes  and  1.8  m  auxiliary  telescopes.  The 
upcoming  near-infrared  VLTI  instrument  AMBER  will  allow 
interferometric  measurements  in  the  near-infrared  J ,  H ,  and 
K  bands  with  a  spectral  resolution  of  up  to  10,000  and  will 
provide  closure  phases.  This  will  enable  us  to  directly  measure 
the  continuum  stellar  diameter  and  to  study  the  conditions  on 
the  stellar  surface,  including  possible  asymmetries  and  surface 
inhomogeneities. 
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